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Abstract 

We  explore  the  conceptualization  of  biaxial  composite  mediums  through  the  process  of  ho- 
mogenization. Biaxiality  is  found  to  arise  when  the  component  mediums  undergoing  homog- 
enization present  two  non-collinear  distinguished  axes.  Two  possible  sources  of  directionality 
in  the  component  mediums  are  considered:  (a)  topological  and  (b)  electromagnetic.  Exam- 
ples of  these  are  investigated  by  considering  the  homogenization  of  particulate  components 
with  (a)  non-spherical  topologies  and  isotropic  electromagnetic  properties  and  (b)  uniaxial 
electromagnetic  properties  and  spherical  topologies. 

1.  Introduction 

In  the  context  of  electromagnetic  material  properties,  the  concept  of  homogenization  is  both  sci- 
entifically and  technologically  important.  Composite  mediums  with  complex  properties  may  be 
conceptualized  through  homogenizing  relatively  simple  constituent  mediums.  Biaxial  symmetry 
in  such  homogenized  composite  medimns  (HCMs)  is  our  primary  concern  here.  We  build  upon 
the  foundation  laid  by  earlier  studies  of  non-dissipative  dielectric  [1]  and  dissipative  dielectric- 
magnetic  [2]  biaxial  HCMs  (wherein  further  background  details  may  be  foimd)  and  generalize 
to  the  bianisotropic  case.  By  considering  only  component  mediiims  of  the  simplest  forms,  we 
demonstrate  that  an  elaborate  HCM  form  can  arise;  and  through  illustrative  parametric  studies, 
we  delineate  S3unmetries  in  the  HCM  structure. 

2.  Preliminaries 

We  consider  HCMs  derived  from  only  two  (distinct)  component  mediums,  each  being  envisioned 
in  particulate  form;  we  refer  to  them  as  the  host  medium  and  inclusion  medium.  Of  the  many 
formalisms  which  have  been  developed  in  order  to  estimate  the  electromagnetic  constitutive 
properties  of  HCMs,  here  we  adopt  the  Bruggeman  formalism  [3,  4,  5].  The  HCMs  emerging 
from  the  nmnerical  calculations  may  be  characterized  by  the  bianisotropic  constitutive  relations^ 

D(x)  = + (1) 

B(x)  = ^£^^^-E(x) +,i,g^^^.H(x).  (2) 

^Vector  quantities  are  in  boldface  while  dyadics  are  double  underlined.  The  unit  dyadic  is  denoted  by  / and 
(ux,  Uy , u,)  is  the  triad  of  cartesian  unit  vectors.  The  permittivity  and  permeability  of  free  space  are  denoted  by 
eo  and  ^o>  respectively. 
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The  HCM  constitutive  dyadics  have  a biaxial  form  which  we  represent  as  [1,  2,  6] 


%JICM  ”1"  ^rC^mr^nr  "I”  ^nr^mT)  , 


(3) 


(-r  = e,  C,  C, 


where  Cr’*  and  6t’’  are  real-valued  scalars  and  we  describe  the  real-valued  unit  vectors  in  terms 
of  spherical  polar  coordinates  as 


= sin  cos  0^^  + sin  6^^  sin  Uj,  + cos  6^^  , 

(X  = r,z;  K = m,Ti]  r = e,^,C,M). 


(4) 


With  one  exception,  we  consider  here  only  constituent  mediums  with  distinguished  axes  lying 
in  the  xy  plane.  Consequently,  our  calculations  reveal  that  the  HCM  unit  vector  pairs  u^c  and 
always  lie  in  planes  perpendicular  to  the  xy  plane  with  the  xy  plane  bisecting  the  angle 
between  u^c  and  u^c-  The  following  identities  therefore  hold 

^mr  = - ^nr  = 4^r  = (t^T  = 4^^  (x  = r,  i;  T = e,^,C,/x).  (5) 

The  one  exception  occurs  when  we  consider  ellipsoidal  inclusions  of  varying  eccentricity;  in  this 

case  one  distinguished  axis  can  lie  along  the  z axis  and  we  shall  treat  this  as  a special  case  in 

Section  3.1.  Furthermore,  since  all  component  mediums  we  consider  are  reciprocal,  results  for 

the  magnetoelectric  dyadic  C need  not  be  explicitly  presented  as  we  find  f = — C 

=HCM  ^ ^ ^ =HCM  =HCM 

All  graphs  are  presented  with  reference  to  the  key  given  in  Table  1.  A volumetric  proportion 

of  inclusion  medium  to  host  medium  of  0.3  is  taken  for  all  calculations. 


3.  Numerical  Homogenization  Calculations 
3.1  Dielectric  case 

We  homogenize  a host  medium  of  permittivity  = 1.2  J and  spherical  topology  with  an 
inclusion  medium  of  permittivity  = (3  + 3z)  I and  ellipsoidal  geometry  characterized  by  the 
shape  dyadic  = diag(2, 1,7).  For  this  particular  example  (and  no  others)  the  identities  (5) 
do  not  hold;  instead  we  take  O^e  = and  4>me  = </>£’*  and  plot  these  angles  as  a function  of  7 
in  Figure  1.  At  points  where  the  inclusion  shape  becomes  spheroidal,  the  unit  vectors  u^^ne  all 
lie  on  a common  axis  and  the  HCM  becomes  uniaxial.  For  all  values  of  7 we  find  that 
is  diagonal  and  hence  the  HCM  belongs  to  the  biaxial  orthorhombic  class  [7].  This  reflects  the 
fact  that  in  this  case  biaxiality  arises  from  a geometrical  structure  based  on  three  mutually 
perpendicular  principal  axes,  namely  those  of  the  inclusion  ellipsoid. 


3.2  Dielectric— magnetic  case 

Here  we  consider  constituents  in  which  the  distinguished  axes  have  an  electromagnetic,  rather 
than  topological,  origin.  We  homogenize  a host  medium  with  constitutive  dyadics  = 

diag(3, 1, 1)  and  an  inclusion  medimn  specified  by  = (1  + i)  A and  = {2  + i)A  where 


A = 


3 cos^  A + sin^  A 2 sin  A cos  A 0 

2sinAcosA  3sin^A  + cos^A  0 

0 0 1 


(6) 


and  both  host  and  inclusion  mediums  have  a spherical  topology.  Plotted  in  Figures  2 and  3 as 
functions  of  A are,  respectively,  the  defining  angles  for  the  unit  vector  pairs  Ume,ne, mu, nfi  and  the 
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corresponding  permittivity  (og’*  and  6«’‘)  and  permeability  scalars  (oj/*  and  . A uniaxial 
dielectric-magnetic  HCM  results  when  the  distinguished  axes  of  the  constituent  mediums  are 
aligned.  With  the  exception  of  the  special  cases  A = 0,  7t/2  and  tt,  all  eight  angles  0e’ft  and 
and  all  eight  scalars  and  belli  have  distinct  values  and  the  biaxial  HCM  is  of  the 
monoclinic/triclinic  type  as  regards  both  and  [7]- 

We  repeat  the  homogenizations  of  Figures  2 and  3 but  now  with  a more  general  host  medium 
characterized  by  = diag(3+3z5,  l+z(i,  l+i5)  and  with  a fixed  angle  for  the  inclusion 

distinguished  axis  of  A = 50°.  The  angles  are  plotted  against  S in  Figure  4.  Our  findings 
for  this  example  may  be  summarized  by: 


Re  Im  Re  = Pr  Im  r*’ 


JiOSt 


= (t>T  = <t>T^  (T  = e,  m),  (7) 


where  Pe,fi  are  proportionality  scalars.  Thus,  the  biaxial  HCM  structure  becomes  orthorhombic 
with  respect  to  permittivity  (permeability)  when  ratios  of  real  and  imaginary  parts  of  and 
gtnc  ^fj^host  are  equal,  despite  the  distinguished  axes  of  the  constituent  mediiuns  being 

non-perpendicular. 


3.3  Bianisotropic  case 

Finally  we  consider  the  general  bianisotropic  case:  the  homogenization  of  an  inclusion  medium 
characterized  by  = 2(l+i)^,  = 1.5(H-i)  A and  = (l+i)A  where  A = 50°, 

with  a host  medium  described  by  _ (Hag{^+Zi5,  l+z<5,  l+z5);  and 

spherical  topology  is  chosen  for  both  component  mediums.  The  computed  biaxial  bianisotropic 
HCM  structure  is  of  the  generalized  monoclinic/triclinic  type  with  angles  0^’^^  and  and 

scalars  and  all  taking  distinct  values,  in  general.  The  corresponding  polar  HCM 

unit  vector  angles  are  displayed  as  functions  of  S in  Figure  5 (the  azimuthal  angles  behave 
similarly).  At  the  point  <5  = 1 we  find  that  the  HCM  is  orthorhombic  biaxial  with  respect  to  all 
four  constitutive  dyadics  and,  as  in  Section  3.2,  the  orthorhombic  state  is  not  associated  with 
perpendicularity  of  the  distinguished  axes  in  the  component  mediums. 
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Tablel:  Key  for  Figures  1-5.  a = a,  5,0  or  (f>. 


Figure  1:  HCM  angles  0^**  and  vs.  inclusion 
ellipsoid  semi-axis  7.  { X = 0^  <f>in  Table  1). 


Figure  2:  HCM  angles  0^*^  and  vs.  inclusion  Figure  3:  HCM  scalars  and  vs.  inclusion 
distinguished  axis  angle  A.  { X = 0^  <f>  in  Table  1).  distinguished  axis  angle  A.  ( = a,  5 in  Table  1). 


Figure  4:  HCM  angles  0g>*  vs.  host  medium  pa-  Figure  5:  HCM  angles  vs.  host  medium  pa 
rameter  5.  ( X = 0 in  Table  1).  rameter  <5.  ( X = 0 in  Table  1). 


